Mechanical properties of the radio frequency field emitted by an antenna array 
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I. INTRODUCTION 

As Maxwell showed, any interaction between electromag- 
netic (em) radiation and matter is inevitably accompanied by 
an exchange of momentum. In addition to linear momentum, 
Poynting |l|] inferred from a mechanical analogy that circu- 
larly polarized light must carry spin angular momentum and 
Allen et al. JUl argued that Laguerre-Gaussian laser modes 
must carry orbital angular momentum. The exchange of spin 
angular momentum between radiation and a mechanical body 
was demonstrated experimentally by Beth in 1935 [3], while 
the transfer of em orbital angular momentum to mechanical 
orbital angular momentum of particles was experimentally 
demonstrated only recently flylla]- 

During the past decade the application of em spin and or- 
bital angular momentum has come to the fore in optics J3l 
and in atomic and molecular physics [18!]. The orbital angular 
momentum of electromagnetic beams now takes an important 
position in various fields of research and applications, ranging 
from manipulating and orien ting small particles [9], astron- 
omy and astrophysics [ 1 Ol 1 1 lL 1 1 2L 1 1 3fl . novel wireless commu- 
nication concepts lfl4lll5h . quantum entanglement [16], cryp- 
tography, and quantum computation [17]. Applications for the 
use of the orbital angular momentum in the radio domain, al- 
lowing digital manipulation under software control, were pro- 
posed in Re f.ll 8L in clu ding radio astronomy and space physics 

Ei m eb eis m nf m, m m mi . 

In the literature, different views can be found concerning 
the possibility of separ ating the spi n and orbital angular mo- 
mentum; see Refs. [3CJ [3j], [33, L3J. In computing the elec- 
tromagnetic beams emitted by a circular antenna array, and 
expressing the angular momentum density explicitly in terms 
of helicity and vorticity, the present paper shows that the an- 
gular momentum of the emitted beams can be clearly sepa- 
rated into parts that belong to helicity and vorticity, respec- 
tively. The separation is associated with the separation of an- 
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gular momentum into spin and orbital parts, a result that has 
been previously shown for optical beams OOl l34tl . Crucial 
for the analytic computation is the fact that the radiated densi- 
ties of energy, linear momentum, and angular momentum are 
homogeneous functions of distance. This allows to neglect 
certain terms in the expressions for these quantities. Since 
helicity and vorticity are uniquely linked to the phases of the 
transverse field components, it is possible to emit and detect 
em beams with spin and orbital angular momentum via trans- 
verse field triangulation, dealing solely with first order quanti- 
ties which can easily be utilized without the need to measure 
weak longitudinal components. 

In order to simplify the present analysis, an approximation 
is made in the geometrical setup by assuming an array with a 
continuum of infinitely many elemental antennas l35ll . 



II. SETUP 

The vector potential of a single or crossed electric dipole 
antenna that is in harmonic oscillation with frequency a> can 
be expressed as 



A<r,)~gj«£ 



(1) 



cf. Eqs. (9.16) and (9.27) in II36I1 . Here i is the imaginary unit, 
k is the magnitude of the wave vector, f is the unit vector in 
the direction of r, Hq is the magnetic permeability, d is the 
length of each dipole antenna, and J is the antenna current. 

Taking the potentials and fields to be complex valued, the 
harmonic time dependence is attributed by multiplying a sta- 
tionary antenna current I with the appropriate complex factor, 



J(f ) = le 



(2) 



The physical fields result from taking the real parts. 

Placing /V crossed dipole antennas on a circle at the posi- 
tions 
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FIG. 1: Plots of the N crossed dipole antennas and their currents. 
The left panel displays the positions and orientations of the crossed 
dipole antennas, and the right panel displays the current vector at 
each crossed dipole for N = 16, h = 1, / = —2, and t = 0. 

where R is the radius of the antenna array [see Fig. Q] (left)], 
and feeding the nth antenna element with the current 
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n=l...N, 



(4) 



[see Fig. [T] (right)], the total vector potential reads 




with 



and 



a/=- 



ikr 
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(5) 



(6) 
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where / is the amplitude of the antenna currents, h = ±1 is 
the helicity associated with circular polarization, and / G Z is 
the topological charge. The topological charge attributes an 
azimuthal phase shift of exp(2niln/N) to the antenna currents 
and is associated with vorticity. 

If the number of antennas N is sufficiently large (typically 
N > 16), the sum in ^ can be approximated by a Riemann 
integral and expressed by its residue in the complex £ plane, 
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In '€ = real 



(8) 
(9) 

(10) 
(11) 



which implies that c ^* c ^ is independent of (j). 



III. FIELDS 

From the vector potential follows the magnetic field 

B = VxA. (12) 

Separating the magnetic field into its far and its near field 
terms, 



B Bftr + B 



far ~r J>near > 



results in 



and 



Bfar = ikr x A 
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(13) 



(14) 



(15) 



where e^' kr Bf ar and e~' kr B neal are homogeneous functions in 
r of degree — 1 and —2, respectively. 

The electric field outside the sources is given by 



E = -VxB, 

k 



(16) 



where c is the speed of light. The electric field results in 

E = Ef ar + Ei nte r + E ne ar (17) 

with 

E far = /fcc[A-(f-A)f], (18) 



Eh 
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and 



Enear = t[V(V ■ (Ae"*-)) - V 2 ■(Ae"*')]^ 



(20) 



The individual terms of the electric field are called the far, Ef ar , 
the intermediate, E; nter , and the near electric field, E nea r. Apart 
from the complex phase, they are homogeneous functions in r 
of degree —1, —2, and —3, respectively. 



IV. DENSITIES OF ENERGY AND MOMENTA 

Imagine a set of concentric spheres centered on the radia- 
tion source at the origin. If a physical quantity Q is conserved, 
the total flow of Q through the sphere surfaces must be the 
same for each sphere. This is only possible if the r 2 depen- 
dence of the surface area cancels with an r~ 2 dependence of 
the density of the conserved quantity. Hence, the density of 
any conserved quantity Q that propagates radially with a con- 
stant speed is required to be homogeneous in r of degree —2. 

The em field energy density is defined by 



u=j(E-E*+c 2 B-B*), 



(21) 
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where Co is the electric permittivity of free space. The conser- 
vation law of energy imposes the condition that the radiated 
energy density must be homogeneous in r of degree —2. This 
condition is fulfilled by the leading order term of the energy 
density, 



"far = ~£ (Efar ■ E^, + C 2 Bf ar • B^) . 



(22) 



The linear momentum density is given by the Poynting vec- 




tor 



FIG. 2: Plots of the function Notice the influence of I on the 

graph. Here Z = (left) and / = ±1 (right). 



p=^9t[ExB*], 



(23) 



which can be split into homogeneous functions in r of degree 
—2, —3, and higher order terms, 



P — Pfar + Pnext + Prest 
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where f , 0, and <j> are the unit vectors of spherical coordinates. 
Since the density of angular momentum scales with distance 
like r~ 2 , the angular momentum is transported all the way out 
to infinity, a fact that is not always well appreciated. 



and 



V. SPIN AND ORBITAL ANGULAR MOMENTUM 



p rest = 0{l/r 4 ). 



(27) 



Due to the conservation of linear momentum, the radiated lin- 
ear momentum density is given by the far field terms, Ef ar and 
B far , only. 

The angular momentum density is defined by 



r x p. 



(28) 



Because 
radiated 
degree - 
pure far 
mentum 
function 
fact that 
term, 



of the conservation law of angular momentum the 
angular momentum density is homogeneous in r of 
-2. This has important consequences. Namely, the 
field does not contribute to the radiated angular mo- 
density, since r x pf ar would then be a homogeneous 
in r of degree — 1, but this term vanishes due to the 
Pf ar is directed along r. Instead, it is the next order 



Jfai 
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Separating the angular momentum density in terms that de- 
pend on the topological charge I only and in terms that depend 
on the helicity h, 



J far Ifar ~H Sfai 



(33) 



allows to identify the topological charge (vorticity) with the 
orbital angular momentum density 

^ ( ^ ( l+c^ 
4a>jU v 8%r ) \ sin0 / 

and the helicity (polarization) with the spin orbital angular mo- 
mentum density 

Sfar = ( 2 (l Sin d(irif) ~ COS e^(*r # ) 

(35) 



that yields the radiated angular momentum density. 

Plugging the vector potential expression © into the expres- 
sions for the magnetic and electric fields and evaluating the 
radiated densities of energy, linear momentum, and angular 
momentum results in 
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Surprisingly, the orbital angular momentum density is eas- 
ier to describe than the density of the spin angular momentum. 
The density of the orbital angular momentum is completely 
independent of the helicity, but the spin angular momentum 
density depends via 'tf*^ on the topological charge. If the 
topological charge vanishes, c €* c € has its maximum at 9 = 0, 
while for any non-vanishing topological charge, c €* c € is zero 
at 0=0, 
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FIG. 3: Radiation patterns of the linear momentum for radio beams 
generated by a circular antenna array with radius 2X; all antennas 
are directly over ground. Without ground the electromagnetic waves 
would also be emitted symmetrically into the lower half space. The 
patterns are for / = (left), / = ±1 (middle), and / = ±2 (right). 




FIG. 5: Radiation patterns for the spin angular momentum. The di- 
rections of the spin angular momentum density vectors alternate be- 
tween — Qh and +9h as displayed by the arrows I* and *1. The patterns 
are for positive helicity, h = +1, and topological charge 1 = (left), 
I = ±1 (middle), and / = ±2 (right). 




FIG. 4: Radiation patterns for the orbital angular momentum. The 
direction of the orbital angular momentum density vectors is in the 
direction of —8/ as displayed by the arrows [* and *]. The patterns are 
for / = (left), / = 1 (middle), and / = 2 (right). 



sign from cone to cone. Hence, the directions of the spin an- 
gular momentum density vectors alternate between — Qh and 
+Qh. 

Often, spin is identified with polarization. We emphasize 
that this is not allowed for the spin angular momentum density; 
only the total spin angular momentum can be identified with 
polarization. For instance, taking I = and h = 1, the spin 
angular momentum density vanishes on the z axis, see Fig. [5] 
(left), whereas the electric and magnetic fields are maximal on 
the z axis and are right-hand circular polarized. 



and 



VII. FLUXES 



V)| 9=0 = 0. 



(37) 



Figure |2] shows graphs of For all plots, an array radius 

of R = 2X was used, where X =2n/k is the wavelength. 



VI. RADIATION PATTERNS 

Radiation patterns of the energy (linear momentum) are 
shown in Fig. [3] If the topological charge vanishes, i.e. I = 0, 
the particular antenna array chosen here emits a beam that has 
maximum intensity on the z axis. If the topological charge 
is non-vanishing, / ^ 0, the fields of the beam cancel on 
the z axis, resulting in a doughnut shaped beam profile that 
becomes wider for larger absolute values of the topological 
charge. 

The far-field approximation of the Poynting vector is pro- 
portional to the energy density and points from the antenna 
array radially outward, 
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(38) 



Therefore, the radiation pattern of the linear momentum is 
identical to that of the energy. 

Radiation patterns of the orbital and spin angular momen- 
tum are displayed in Figs. [4] and [3] respectively. The orbital 
angular momentum density always points in the direction of 
—0/, whereas the spin angular momentum density changes its 



The total fluxes of energy, linear momentum, orbital, and 
spin angular momentum follow from integrating the densities. 
Let S be a closed surface with the antenna array inside it. The 
total fluxes through the surface are 



dF 
dt 



fcr 2 dQ. , 



(39) 



where F stands for the total flux W, P, L, or S and / for the 
corresponding density Mf ar , pf ar , lf ar , or Sfar, respectively. Be- 
cause of the rotational symmetry around the z axis, the orbital 
and spin angular momenta, L and S, are both parallel to the z 
axis. The linear momentum is radiated symmetrically in the 
positive and negative z directions. Thus, the net vanishes, 



P = 0, 



(40) 



and the angular momentum of our electromagnetic beams be- 
comes an intrinsic but non-local property. 

Expressing the orbital and spin angular momentum in mul- 
tiples of the total energy, we obtain l2l l30ll 
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and 



W 



dS z /dt 
dW/dt 



1 

CO 



h 

CO 



(41) 



(42) 



5 



Torque 




and 



Absorber 



Ground 

FIG. 6: A radio experiment that measures the torque of the total 
angular momentum flux. 



Jfar 4c0jU V Mr ) \ sin0 / v ; 

It is useful to compare these densities with those given in 
Section HVl Eqs. d30b — (l32t. The essential difference is in the 
angular momentum density. For pure longitudinal currents, 
there is only an orbital part, 



Jfai 

and the spin part vanishes, 



far ' 
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VIII. DETECTING THE TORQUE 

The torque of the angular momentum flux can, at least 
in principle, be directly measured with an absorber that sur- 
rounds the antenna array in the upper half space; see Fig. [6] If 
the antenna array is near ground with reflectivity p g = / re fi /I m , 
and if the absorber has a reflectivity of p a and a transmittance 
of r a = ^trans An, the torque is 



d(L z +S z 
dt 



a- 



1 + hdW 
co dt 



(43) 



where dW /dt is the power that is radiated by the antenna array 
and a is determined by the reflectivities and the transmittance 
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'far 



0. 



(51) 



This shows that electromagnetic waves associated with pure 
longitudinal currents can carry orbital, but no spin angular mo- 
mentum; cf. the discussion of longitudinal plasma modes in 
Ref.lM 

Because of the term (1 — cos 2 6) the beams of electromag- 
netic waves associated with longitudinal currents have a deep 
null along their axis of propagation which gives rise to a 
doughnut shaped beam profile, even in the absence of any 
topological charge. 

The LOIS radio facility [HHHD], currently under construc- 
tion in Sweden, is the only device that utilizes 3D vector sens- 
ing antennas and is readily able to detect these waves. Other 
radio telescopes, e.g., LOFAR [40] and SKA [41], can be eas- 
ily reconfigured to achieve the same capability. 



RADIO WAVES ASSOCIATED WITH LONGITUDINAL 
CURRENTS 



IX. 



Based on an idea put forward by Carozzi [37], let us re- 
place the transverse crossed dipoles of the circular antenna ar- 
ray by electric dipoles that point in the longitudinal direction 
and feed the longitudinal dipole antennas with the currents 




n=\...N. 



(45) 



Then the vector potential reads 




A = — Mr)Sf(0,0), 



(46) 



where srf and ^ are as given in Sectionllll 

Evaluating the densities of energy, linear momentum, and 
angular momentum, we obtain 
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X. GENERALIZING THE RESULTS 

From the previous Sections one might get the impression 
that the analysis is specific to radio beams emitted by a circu- 
lar antenna array. However, as we now will show, the analysis 
is more general and device independent. 

Any device or source that generates a current distribution 
on a circle 



where p, < t > , and z are cylindrical coordinates, results in the 
vector potential 
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cf. Eqs. (O-®. Hence, for the current distribution given here, 
the results presented in the previous Sections are valid. 

Since the angular dependence of the vector potential can be 
expressed by its residue, Eq. (O, there is no need for the cur- 
rent distribution on the circle itself. Mathematically speaking, 
the current distribution can be along any path that encircles the 
origin, provided that the residue is kept unchanged and that 
no further vortices contribute. Moreover, the current distribu- 
tion can be chosen arbitrarily in the xy plane, provided that 
the residue is kept unchanged. Even if the residue changes, 
assuming that Eqs. ( fTOb and dTTb remain valid, all results pre- 
sented in this paper, except Figs.[2]-|5] will hold. 

If the current distribution is not in a plane, but in three- 
dimensional space, knotted optical vortices appear [42] . 



XI. INTERFERENCES BETWEEN WAVES ASSOCIATED 
WITH TRANSVERSE AND LONGITUDINAL CURRENTS 



respectively. 

If integrated, most of the interferences vanish. This allows 
to carry over the principle of superposition of electromagnetic 
waves to their mechanical properties, 
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(65) 
(66) 
(67) 
(68) 



However, for certain quantities the interferences do not can- 
cel: 

4 nter V0, Lj ntel V0, s» terf ^0, S; nterf ^0. (69) 
The consequences of this is the subject of a forthcoming paper. 



If the vector potential reads 

1 /«\ 
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the densities of energy and momenta are 
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The contributions associated with transverse currents are here 
denoted with the upper index "_L" and are given in Sections 
[lV]and[V] Eqs. (f30b — (f35b. The contributions associated with 
longitudinal currents are given in Section IIXI and the contri- 
butions due to interferences read 

" interf - T7T ( TTr) 2 sin 9 cos cos (A0 + arg * a) , 

(61) 
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XII. DISCRETE SOURCES AND THE UNCERTAINTY 
PRINCIPLE 

If, as in Section [Xj the sources are continuous, the analy- 
sis is exact in the far field and the number of orbital angular 
momentum states is unbounded. In contrast, if the sources are 
discrete, the presented analysis is an approximation. Luckily, 
the approximation improves very fast with the number of dis- 
crete sources, and there is typically no need to worry about 
the accuracy. 

However, for the case that only very few discrete sources 
are present, or if one is interested in the number of orbital an- 
gular momentum states that can be utilized, we propose an al- 
ternative approach. Namely, starting from a continuous source 
distribution the discrete nature of the sources is imposed via 
a "Pacman" shaped mask that has slits at the appropriate po- 
sitions. Due to Fourier convolution with the mask, the un- 
certainty principle enters and gives rise to a distribution in the 
orbital angular momentum 143114411 . This distribution has a pe- 
riodicity of N and peaks of width proportional to 1 /N limiting 
the number of orbital angular momentum states that can be uti- 
lized in a circular array of N equidistributed discrete sources 
to 



#{1}=N-0(1/N), 



(70) 



where the implied constant in the last term depends on the 
noise level and on the details of the discrete sources. 



XIII. CONCLUSIONS 

Our analysis of the properties of radio beams emitted by a 
circular antenna array highlights analytically profound aspects 
of the beam spin and orbital angular momentum. Namely, how 
the beams are connected to helicity and vorticity; that the an- 
gular momentum of the beams is an intrinsic property; that it 
is transported all the way out to infinity; that it can be sepa- 
rated into spin and orbital angular momentum; and that spin, 
unlike its density, can be identified with polarization. 
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In particular, separating the angular momentum into spin 
and orbital angular momentum, we find that the density of 
orbital angular momentum depends on the vorticity only and 
can be expressed in quite simple terms, whereas the density of 
spin angular momentum depends on both helicity and vortic- 
ity. 

The approach via an antenna array illuminates the physics 
from a new perspective and confirms previous studies of or- 
bital angular momentum of electromagnetic beams j2l [30l[34ll . 
Emerging as a byproduct, but being of interest in its own, is 
the finding in Section llXl 
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